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Abstract 



We consider the biorthogonal polynomials associated to the two-matrix model where the eigenvalue distribution 
has potentials Vi, V2 with arbitrary rational derivative and whose supports are constrained on an arbitrary union 
of intervals (hard-edges). We show that these polynomials satisfy certain recurrence relations with a number 
of terms di depending on the number of hard-edges and on the degree of the rational functions V/. Using 
these relations we derive Christoffel-Darboux identities satisfied by the biorthogonal polynomials: this enables 
us to give explicit formulas for the differential equation satisfied by -|- 1 consecutive polynomials, We also 
define certain integral transforms of the polynomials and use them to formulate a Riemann-Hilbert problem for 
(di -I- 1) X {di + 1) matrices constructed out of the polynomials and these transforms. Moreover we prove that 
the Christoffel-Darboux pairing can be interpreted as a pairing between two dual Riemann-Hilbert problems. 
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1 Introduction and setting 

In this paper we consider the biorthogonal polynomials associated to the two-matrix model. The model is 
defined by a measure on the space of pairs of Hermitean matrices Mi, M2 of the form 

d/i(Afi, M2) := dAfidAf2e-T''(^^(*^^)-^=^(*^^)+*^^*^=) . (1-1) 
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Using Itzykson-Zuber/Harish-Chandra's formula the model can be reduced to the study of biorthogonal poly- 
nomials [14] (BOPs for short), namely two sequences of polynomials {7r„(x)}, {crn(y)} 

f f da;d2/e-^i(^)-^^('')+^^^„(x)a™(y) = 5„„ . (1-2) 

JR JR 

For the model to have a probabilistic interpretation, the potentials should be real and satisfy certain growth 
conditions to ensure the convergence of the integrals. In order to introduce the setting of this paper we 
consider the following situation (which is strictly included in the more general setting to be expounded later) 

1. There is a finite collection of disjoint intervals I = [jlj c M^. and J = [jj Jj c Mj, (Rx denotes the real 
axis of the x-variable), in the complement of which the potentials are -|-oo: in other words the matrices 
Ml, A/2 have spectrum confined to these multi-intervals, so that the associated BOPs satisfy 

da;d2/c-^i(^)-^=('')+^%„(x)(T„(y) = Snm (1-3) 

2. The two potentials Vi{x) and V2{y) are the restriction to I, J (respectively) of real-analytic functions 
with rational derivative (with poles symmetrically placed off the real axis, or on the complement of the 
intervals on the real axis) together with the necessary growth condition if the intervals are unbounded. 

This situation has been addressed in [3] within the general context of bilinear moment functionals. Indeed it 
is convenient to recast the orthogonality condition in a more abstract setting where one considers a bimoment 
functional C : C[x] g)C[?;] C defined by 

Cix'\y^) := J J da;d2/x'y^e-^i(^'-^=(«)+^^ = . (1-4) 

and then extended by linearity to arbitrary polynomials. The biorthogonality condition then reads 

C{TTn\<Jm) = Snm ■ (1-5) 

The properties of the potentials Vi, V2 and the supports of integration can be dealt with on the same footing 
by purely algebraic methods: to this end one introduces four polynomials Ai,Bi, i = 1,2 according to the 
strategy outlined hereafter. Let {xj,mj) be the location of the poles of V-[{x) with their order (we include all 
of the poles, in this case also the complex conjugates, which clearly come in with the same multiplicities) and 
let ttj be the endpoints of /. We define then Ai,Bi (and similar expressions for A2,B2) as follows 

Bi{x) = Ylix - x,T' Y\{x - a,) , := V{B, - B[ , (1-6) 

so that now = — ^g— ^- It is a straightforward exercise to verify (using integration by parts) that the 
bimoment functional satisfies the following distributional identities for arbitrary p{x) e C[x], s{y) e C[y] 

/:( - B,(x)p'(x) + Ai(x)p(x)|s(j/)) = c(^B,(x)p{x)\ys(y)') , (1-7) 
/:(p{x)\- B2{y)s'{y) + A2iy)s{y)) = c(xp{x)\B2{y)s{y)') . (1-8) 
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Abstracting formulas (1-7,1-8) from the specific context, we will say that a bimoment functional L is semi- 
classical if it satisfies those same relations (1-7,1-8) for some given (and fixed) polynomials Ai,Bi. The 
name comes from a similar usage in the context of ordinary orthogonal polynomials [13]. 
Such functionals have been studied in [3], where it was shown that 

Proposition 1.1 For given Ai, Bi, i = 1, 2, a semiclassical moment functional C is the linear combination of 
S1S2 independent functionals C^^^, /i = 1 . . . , si, v = \, . . . , S2, where Si = max(deg A^, degi?i -I- 1) 

More importantly (at least in the case dcgA^ ^ dcgBi -\- 1) all of these moment functionals C^^i, can be 
given an integral representation completely analog to (1-4), but without any restriction on the reality of the 
potentials or of the contours of integration: this is the setting of the present paper. 

1.1 Connection to other orthogonal polynomials 

The algebraic properties of semiclassical bilinear moment functionals apply to a slightly different class of 
orthogonal polynomials. Let us consider in fact orthogonal polynomials in the complex plane with respect to 
a measure of the form 

dM(z,z):=e-l^l'+29?v(.)d2^ (^_9) 

where V{z) is a holomorphic function such that V'{z) is rational. The convergence of the measure mandates 
that the residues of V'{z)dz must have real part greater than — ^ and that the behavior at 00 of y cannot 
exceed the second power (and also a certain open condition on the coefficient of this quadratic term which we 
do not specify here). Orthogonal polynomials are defined as a holomorphic basis of _L^(C,d/i). It is amusing 
to note that the moment functional 

C{z'\z'') := JJz^z'=d/i(2,z) =: ^^-^ (1-10) 

c 

is a semiclassical moment functional (using Stokes' thm. in vece of integration by parts) with just some 
(obvious) reality constraint on the bimoments. Therefore all the algebraic manipulations that rely on the 
semiclassicity alone carry out verbatim to this case and confirm certain manipulations used in [21]. With very 
minor and trivial modifications. Section 2 almost entirely generalizes (in particular Thm. 2.1). Significant 
differences (sufficient to require a different analysis to appear elsewhere) arise in the construction of the 
fundamental systems and the Riemann-Hilbert problem. 

1.2 Connection to 2-Toda equations 

The framework of this paper is connected to the general theory of 2-Toda equations [22, 1, 2]. This is the 
theory of a pair of (semi)-infinite matrices P,Q (in our notation) where Q is lower-Hessenberg and P is 
upper-Hessenberg'^ which evolve under a bi-infinite set of commuting flows {ij,ij}jgN 

dt,Q = -j[Q , (Q^)-o] , %Q = - j[Q, {P')-o] (1-11) 

^We say that a matrix is lower Hessenberg its (i, i + 1 + k) entries vanish (Vfe = 1,2,...) and also all (i, i + l)-entries are 
nonzero. A matrix is upper Hessenberg if its transposed is lower-Hessenberg. 
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8t,P = -j[P, {Q')+o] , d^ P = --[P, (P^)+o] (1-12) 

where the subscript +0 denotes the upper/lower triangular part plus half of the diagonal (we are assuming the 
normalization such that the upper triangular part of Q coincides with the transposed of the lower-triangular 
part of P). 

Let now Q,P be semi-infinite matrices. We can use Q,P to denote the matrices expressing the multi- 
plicative recurrence relations of a sequence of polynomials, 

n+l n+1 
XTTn = 2 QnjTTj , 2/CT„ = ^ P„jCr„ , (1-13) 

where the polynomials are recursively defined by this relation. Using the generalization of Favard's theorem 
proved in our [3] we prove the existence of (unique) a bimoment functional C : C[x] ®C[?;] C such that 

C{TTn\arn) = <5nm • (1-14) 

It then follows easily that the 2-Toda flows are linearized by this moment map, in the sense that the 
solutions (3(t, t), P(t, t) are simply the multiplication matrices for the biorthogonal polynomials of the moment 
functional 

C,rt(.»\») ■■= ^ . \e-^^y\'^ . (1-15) 

The moment functionals of semiclassical type (eqs. 1-7, 1-8) that we are going to analyze form a particular 
class of reductions of the above-mentioned 2-Toda hierarchy. The simplest situation is the one of bimoment 
semiclassical functionals with polynomial potentials as the ones considered in [4, 5, 6], where the matrices 
P, Q are also finite band. Moreover the solutions which arise in the context of semiclassical bilinear functionals 
also satisfy the (compatible) constraint of the string equation 

[P, Q] = hi (1-16) 

(the constant h can be disposed of by a rescaling). The parameters of the (finite-dimensional) reduction are 
the coefficients of the potentials: for more general semiclassical moment functionals as the ones considered in 
this paper, the parameters involve not only the coefficients of the partial fraction expansions of the (derivatives 
of the) potentials, but also the position of the poles and the position of the end-points of the supports of the 
measure (the hard-edge endpoints). 
The paper is organized as follows 

1. In Section 2 we derive the recurrence relation satisfied by the biorthogonal polynomials of a semiclassical 
moment functional. There are two types of recurrence relations: one which involves the multiplication 
by the spectral parameter (and plays the role of the more standard three-term recurrence relation for 
orthogonal polynomials) and one which involves a differential operator acting on the polynomials. 

2. In Section 3 we recall some possibly not well known facts about a certain class of linear homogeneous 
ODEs. These equations are next in simplicity to the class of constant coefficients ODEs, inasmuch as the 
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coefficients are allowed to be linear functions of the independent variable. When considering the formal 
adjoint equation then the classical bilinear concomitant provides a nondegenerate pairing between the 
solution spaces of the pair of mutually adjoint ODEs. In this case we give an interpretation of it in terms 
of an intersection pairing between certain contours used in the representation of the solutions as 
contour-integrals. This part of the paper is logically quite independent on the rest but it is nevertheless 
necessary in order to understand certain constructs of the following section. 

3. In Section 4 we define the auxiliary wave vectors for our functionals, using a certain multiple integral 
transform which relies upon the form of the bilinear concomitant associated to our semiclassical moment 
functional (extending some of the results of [8]). These expression will prove crucial in the formulation 
of a first order ODE of rank di + 1 = 1 + dcg{Ai) satisfied by the biorthogonal polynomials. We also 
derive the analog of the Christoffel-Darboux identities satisfied by standard orthogonal polynomials to 
our case of biorthogonal polynomials: similar expression were extensively used in [5, 8] for the case 
where the potentials Vi are polynomials (which is a subcase strictly included in our present setting) and 
in absence of hard-edge endpoints. The novel feature is that these new identities involve not only the 
biorthogonal polynomials of the moment functional C itself, but also those of the associated bilinear 
semiclassical moment functionals 

£:=£(Bi.|.); £:=/:(.|B2«) • (1-17) 

This feature appears prominently in the perfect duality of the Riemann-Hilbert problems appearing in 
the next section. 

4. In Section 5 we define a pair^ of piecewise-ana lytic matrices constructed out of the entries of the wave- 
vectors and their auxiliary wave-vectors. They satisfy certain jump conditions on contours in the complex 
plane and some asymptotic behavior at the zeroes of Bi. Moreover they satisfy rational first order ODEs 
with poles at the zeroes of Bi. The Christoffel-Darboux identity, when written as a bilinear expression 
for these matrices becomes a perfect pairing (Thm. 5.1) in the sense that establishes a nondegenerate 
constant (in x) duality-pairing between the two solution spaces. This pairing should be thought of as 
the "dressed" form of the bilinear concomitant pairing introduced in Sect. 3. Similar Riemann-Hilbert 
problems have appeared elsewhere in the literature, e.g. [19, 18, 8, 4]. 

In order to facilitate the navigation through the paper all proofs of more technical nature are collected in the 
appendix and only those that may help the understanding are left in the main body of the paper. 

Acknowledgments. The author would like to thank R. Teodorescu for discussion during the summer 
2005 conference on Random Matrices at CRM, John Harnad for daily stimulating interaction and one of the 
referees for pointing out some mistakes in an early version. 

*In fact there are two such pairs, the other being obtained by interchanging the roles of x,y, Bi,B2 etc. 
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2 Semiclassical bilinear moment functionals of type BB 



We consider an arbitrary bilinear semiclassical moment functional (as defined in the introduction) [3], i.e. 
satisfying (1-7, 1-8). Let qi = deg{Bi) and dj = deg(Ai): we assume that di ^ qi + 1 ("type BB" in the 
terminology of [3]). We also make the assumption that the two pairs of polynomials Ai,Bi are reduced 
in the sense that the only common zeroes of Ai and Bi [i = 1,2) are amongst the simple zeroes of Bi. 
Any moment functional coming from a representation like the one in the introduction (1-4) has this property 
of reducedness. In [3] the case of non-reduced moment functional is also considered, and it corresponds to 
functionals which may be expressed as delta functions (or derivatives thereof): we refer ibidem for details. 
It is known [3] that any such reduced moment functional can be expressed in integral form 

di d2 

My := Cix^y^) =2 I] ^^..^^^M^'ly') (2-1) 
C^uix'W) = f r e-^i(")-^=(«)+"Ma;dy (2-2) 

Vl(y) = (2-3) 

£ = JJ x'y^c-^^'^'=^-^^'^y^+'=ydxdy (2-4) 

(2-5) 

The two sets of contours of integration r^^,^ and Ty i^ are defined in the x and y complex planes respectively 
and in completely parallel fashion: we will define them in Section 3.1. We have also introduced the short-hand 
notation 

Note that the case of hard-edges is included: the hard-edges are the zeroes of Bi that cancel with the zeroes 
of the denominator defining V{ in eq. (2-3). 

The constants k^,i, e C are arbitrary (not all zero). In the paper we will often invoke "genericity" conditions 
for the moment functional C: by this we mean that the genericity is in the choice of the >f-constants and not 
in the choice of Ai, Bi which we consider as given once and for all. All of the genericity conditions that we will 
use can be translated into the nonvanishing of certain infinite sequences of minors of the matrix of bimoments 
M = \pLij\- since the moments are linear in k as per (2-1), this genericity boils down to avoiding an 
at-most-denumerable collection of divisors of homogeneous polynomials in the >f-space. 

2.1 Biorthogonal polynomials 

Let us consider the biorthogonal polynomials associated to this bilinear moment functional, namely two se- 
quences of monic polynomials satisfying the following conditions 

{7r„(a;),cr„(j/)}„eN 
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7r„(x) =x" +0(x"-i) 
a„(y)=2;" + 0(2/"-i) 

'C(7r„|(T,„) = hnSnrn ■ (2-7) 

The existence of these BOPs is guaranteed provided that the principal minors of the matrix of bimoments do 
not vanish 

A^jC] := det[My]o«;^j,^:„-i ^0 Vn e N , (2-8) 

which also guarantees that hn # , Vn e N ([3]). We find it more convenient to deal with the normalized 
BOPs; 

Pn ■■= -j= , Sn := — ^ (2-9) 
We will use the following quasipolynomials 

:= p„e-^^(^) , </.„ := s„e-^^(^) (2-10) 

and the following semi-infinite vectors (wave vectors) 

p{x) ■.= [po,pi,...,p„,...Y , s{y):=[so,si,...,Sn,...f (2-11) 
^:=p(x)e-^i(") , $ := s(?/)c-^^(«) (2-12) 

It will become necessary to consider the following associated semiclassical functionals defined by the relations 

Cip\s) := C{p\B2 s) , C(p\s) := £(Bi p\s) . (2-13) 

We leave to the reader the simple check that these are also semiclassical moment functionals where the 
potentials are replaced -respectively- by 

I V2{y) ■.= V2{y)-lnB2{y) 
\v2{y) = V2{y) 

These definitions amount to A2 = A2 - B'^, B2 = B2 so that = ^^gf" = ^■ 

Note, however, that they are defined along the same contours as C and with the same coefPicients >c's. 

2.2 Multiplicative recurrence relations 

We now prove 
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Theorem 2.1 The BOPs satisfy the following finite-term recurrence relations: 



92 



)d2 
J=-l 

/ 91 \ di 

qi = deg{Bi), d, = deg(Aj) , 



(2-16) 
(2-17) 



where ij{n) = for n ^ d2 and mj{n) = for n ^ di, under the genericity assumption (to be further 
discussed In Remark 2.1) 



A„ 9 := det 



MlO • • ■ Ml. 92-1 
M20 • • ■ ^2,92-1 



^",92 — 1 



Moo 
Mio 



MO.n — 52 — 1 
Ml.n — 92 — 1 



^ , Vn > 52 6 N. 



(2-18) 



M^^ — 1,0 ■ ■ ■ M^ — 1,?^ — 92 — 1 

T/)e coefficients q;_i (?i) anJ (n) are nonzero for any n; furthermore, under the same genericity assumptions 
letting Gi, hi be the leading coefficients of Ai, Bi we have 

b2ad2 (n)V^„-d2 = a2(q2 {n)^Jhn-q^ ^ , d2 

biPdt {n)^/KZZ = 01^91 {n)^Jhn-q^ Q , n ^ di (2-19) 

Proof. We prove only one relation, the other being proved by interchanging the roles. 

The statement a_i{n) ^ follows from the form of the recurrence relation by comparison of the leading 
coefficients, which gives 

a_,{n) = .^^0. (2-20) 



The fact that (j{n) = for n ^ c?2 is a choice of convenience: indeed, since d2 > 92 any xp„ can be written 
as a linear combination of the same BOPs of degrees m = 0, . . . ,n + 1 for n ^ d2. 

Consider xpn{x): by "integration by parts" (i.e. using relation 1-8 from right to left), we immediately conclude 
that 

xpnix) 1 i32(y)C{l,y,...,2/"-^=-i}=:Fi2) (2-21) 

Therefore v}^'^^^ is in the common annihilator of xpn{x), . . . ,xpn-q2{x). We now show that it is generically 
possible to fix the coefficients £j{n) of a linear combination as the left hand side of eq. (2-16) such that the 
result is perpendicular to any polynomial q{y) of degree dcg(g) < n — d2. Let 



q{y) =B2{y)a{y)+b{y) 
be the long division of q by B2 with remainder b: then 

C(xpn{x) q{y) ) = C[xpn{x) B2{v)a{y) +b{y) ] = jO[ xpn{x) b{y) 



(2-22) 



(2-23) 
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Since the remainder b{y) is of degree at most q2 — 1, we can find the aforementioned linear combination by 
solving the system 



y 



fc = 0, . . . , 52 - 1 



(2-24) 



After doing so we have that a suitable linear combination in xC{pn, ■ ■ ■ ,Pn-qi} is perpendicular to any 
q = + b with dcg(a) < n — d2 — 92. dcg(5) ^92 — 1, or -in other words - to any q{y) of degree less than 
n — c?2, thus proving the shape of the recurrence relation. 

In order to clarify the genericity assumption we are imposing we express the above condition as a non- 
vanishing condition of certain submatrices of the matrix of moments. Indeed the polynomials p„ := p„ -|- 
^j{^)Pn-j SI'S uniquely determined by the condition that (for ?? > (72) 



1. The degree of is n; 

2. The polynomial p„ is ^-orthogonal to 1,?/, . . 

3. The polynomial xpn is ^-orthogonal to l,y, . 



■ ,y 



1/92-1 



This determines them (for n > (72) as the following determinants (up to a nonzero multiplicative constant) 



Pn 



det 



Mio 
M20 



-1-1,0 



^1,92-1 
^2,92-1 



A'n-H,ij2-1 Mn.O 



Moo 
Mio 



A^O.n — (72 —1 
Ml.n — g2 — 1 



1 

X 



1 X" 



(2-25) 



The genericity condition is then the nonvanishing of the principal minor of size n of the above expression, 
namely the nonvanishing of the determinants advocated in the statement of the theorem (eq. 2-18). 
The normalization that Pn = Pn + {lower degree) gives for the c„ of eq. (2-25) 



1 



A„ 2 



(2-26) 



Let us now check that this genericity assumption is actually equivalent to requiring (n) ¥= , Vn. Denoting 
by 02,^12 the leading coefficients of A2{y), B2{y) we find 



b2adAn)VhZ:Z = C{xpn\B2y''-^'-^-) = C(pn\A2y''-'''-'^' -Oiy^-"--^)) = £(p„|a2y"-«^) = e,,in)a2^. 

(2-27) 
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This proves the identity (2-19): to prove that it does not vanish under our genericity conditions we compute 



/:(p„|a22/"-'^) 



a2 



:det 



Mio 
M20 



Mn+1,0 



Moo 

MlO 



Mn,0 



MO,n — g2 
Ml,n — g2 



IJ'n,n — q2 



02 A„ 



A„.2 vft^ 



^0 



(2-28) 



Q.E.D. 

We can represent the previous recurrence relations in matrix form as follows 
Proposition 2.1 T/je wave vectors satisfy the following recurrence relations 



(2-29) 



where L is the lower triangular matrix with (72 subdiagonals whose matrix entries are Lnm = tn{n — m) and 
A Is a lower Hessenberg matrix with entries Anm = Q;„(m — n) (similarly for M, B ). The entries in the lowest 
and highest diagonals in 1 + L,A are non vanishing. 

2.3 Differential recurrence relations 

Proposition 2.2 Under the genericity assumption^ that the principal minors of the associated moment func- 
tionals C, C are all non-vanishing (or -which is the same- the existence of biorthogonal polynomials for C, C), 
the BOPs satisfy the following differential finite-term recurrence relations 

/ qi \ '''^ ^ 



Vy ^S„ -I- ^ ijin + j)Sn 

V, := 8^ - V{ix) , 



= 8y - V^iy) . 



(2-31) 
(2-32) 



In matrix form we have 



V,(l +M*)p = -B*p 



where the matrices above are defined by 



Mnk = ffln-k{n) , Bnk = (3n-k{n) 

Lnk = in-k{n) , Ank = a„_fe(n) . 



(2-33) 



(2-34) 



Note that they have the same shape as AI, B, L, A respectively (whence the mnemonics of the symbols). 

5 See Remark 2.1. 
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Proof. We prove only the first of the two relations, the other being proved analogously. Consider the unique 
(generically existing) vector p„ in C[p„, . . . ,p„+gj which is divisible by Bi{x) and "monic" w.r.t. p„ in the 
sense that p„ = p„ + C[p„+i, . . . ,Pn+qi]- Writing then p„ = BiQn we find 



e^'d^Pne-^' = B[qn + B[q'^ - V{Biqn = Biq',^ - Aiq„ 



(2-35) 



This implies that {—8x + V()p„ is a polynomial of degree n + di in spite of the fact that ¥( is rational. 
Moreover 



C[{-dx+Vi)pr. 



V 



d -Biq'^ + Aiq,, 



y 



jOi Biqn 



k + l 



k + l 







' J 

k<n-l (2-36) 



Note that the above relation (2-36) is implicitly an assumption on the existence of polynomials qn of exact 
degree n which are £-orthogonal to all lower powers of y: this is equivalent to saying that there must exist 
the BOPs for C, whence our genericity assumption in the statement of the theorem. Q.E.D. 

For later convenience we also remark that the genericity condition we are invoking now is also equivalent 
to requiring that the vectors (the superscript denoting the r-th derivative) 



Bi{xj) = 0, r = 0. . .rj , Bi{x) = h ]^(a; - Xjf 



(2-37) 
(2-38) 



be linearly independent: indeed 



91 

Pn + 2 fhj (n + j)Pn+] = er, 



Pn{xi) 



Pn + qA^l) 



Pn''\xs) 



Pn+qi i^s) 



Pn{x) Pn + l{x) 



Pn+qi {x) 



(2-39) 



where e„ is the inverse of the {qi + 1, l)-cofactor of the above matrix. The proposition can be rewritten for 
the wave vectors as follows 



Proposition 2.3 The wave vectors satisfy the following differential equations 



(2-40) 



where Mnk = frik-nin) and Ank = oik-n{n) (and similar expressions for AI, B). 

The matrices M, B, L, A play the same role of M, B and L, A for the moment functionals C and C respectively. 
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Proposition 2.4 The vectors of polynomials ' 



p(x) := (l+L)-ip , %) 



B2iy) 



(1+L*)s(y) 



(2-41) 



(where L (and M) are defined by eqs.(2-33) of Prop. 2.2) are the biorthogonal polynomials for C. Similarly 
the vectors of polynomials 

1 



p{x) 



Bi{x) 

are the biorthogonal polynomials for C 



(1 + M*)p, s(y) := (l + A/)-is 



(2-42) 



Proof. The two statements are completely parallel and hence we prove only the first. 

By definition of the matrix L in Prop. 2.2 the polynomial entries of (1 -I- L*)s are all divisible by B2 
therefore s is indeed a vector of polynomials. Next we have (using an obvious matrix notation) 



s* (1 -fL) = 1 Q.E.D 



C (p s*^ =c(^{l + L)-^p s*(l + L)^ = (1 + L)-i£ (p 
We also have 

Lemma 2.1 The matrices L, A, AI, B and the matrices L, A, M, B are related by 

A{1+L) = {1 + L)A , B{l + M) = {1 + M)B . 
Proof. Once more we prove only the first. 



(2-43) 



(2-44) 



C(xil+L)p 



A{l+L) = C( Ap 



= £ (1 + L)p 



s\l+L) = £ x{l+L)p 



{-B2dy+A2)s*] =£ (1 + L)p 



£ (1 + L)p 



Vys\l + L) =Ci{l+L)p 



s\l + L) 

- VyB2s' 

s^a] = (1 +L)A Q.E.D. 



(2-45) 



Lemma 2.2 The associated wave vectors p,s and p, s satisfy 

x{l + L)p = Ap 

y{l + M)s = Bs (2-46) 
Moreover, under the same genericity assumptions 

m,,(7i) #0#4(n) Vn (2-47) 



h2ad^{n)\l hn-d2 = 02^92 ('^)'\/ ^"-92 ^ ^ (2-48) 



®The expressions {\ + L) ^ etc. are defined by the geometric series; since L is strictly lower triangular, such geometric 
series is entry-wise well defined. 
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Functional 




BUFs 


TV /T li- 

Mult, rcc 




Din. Kec. 




£(.|.) 




p,s 


x{l +L)p = 
y(l + M)s = 


Ap 

Bs 


V,(l + M*)p = 

Vy(i + £*)s = 


-S*p 

-l*s 


C{.\.) = C{.\B2») 


p ■ 

s :-- 


= (1 +Z)-ip 


x{l + L)p = 

■k 


Ap 


Vj.(l + M*)P = 
* 




£(.|.) = /:(i?i.|.) 


P = 

s = 


^{1+M*)p 
= (1 + M)-is 


* 

y{l + M)s = 


Bs 


Vy(l+Z*)s = 


-I*s 






^(1 + M*)p 










£(.|») = /:(Bi • 1^2.) 


s = - 


(l + £)-^P 
\ (l + M)-is 


x{l +L)p = 
y{l + M)s = 


Ap 

Bs 


* 





(1 + = ^(1 + £) (1 + £)! = 1(1 + £) (2-50) 

(1 + M)B = B{1 +M) (1 + M)B = B{1 + M) (2-51) 



Table 1: the various recurrence relations. 

Proof. Recalling that p = (1 -I- L)^'^p (by definition), we find 

x{l + L)p = xp = {l+ L)-^Ap = A{1 + L)-^p = Ap . (2-49) 

The relations (2-48) for the moment functionals C,C are proved in exactly the same way relations (2-19) are 
proved for C Q.E.D. 

We can summarize all the relations collected so far in Table 1. Here the matrices A, A, A are lower- 
Hessenberg matrices with d2 nontrivial sub-diagonals, B,B,B are lower-Hessenberg with di nontrivial sub- 
diagonals. The matrices L,L,L and M,M,M are strictly lower triangular matrices with q2 or qi nontrivial 
subdiagonals respectively. 

The ★'s mean that there are (possibly under similar genericity requirements for the corresponding functional) 
similar relations as in the corresponding box on the first line, for which however we do not need to define 
symbols for our purposes. 

Remark 2.1 We now address the genericity assumptions invoked in Thm. 2.1 and Prop. 2.2 in the case of 
real potentials and support on the real axes as discussed in the Introduction. For Prop. 2.2 the assumption is 
simply the existence of the BOPs for the associated moment-functional L; in the case of real potentials with 
supports on the real line one can follow [12] and show that BOPs do exist. Since in £ and £ have the same 
supports and B2{y) would be positive on the support (provided that none of the higher-multiplicity zeroes lie 
within the support) then one can conclude that for the case of relevance to the Hermitean two-matrix model 
the Prop. 2.2) is always valid. 

Less transparent is the extent of the limitation imposed by the genericity assumption (2-18) used in Thm. 
2.1; the polynomials pn appearing in the proof of said theorem play the same role in respect to pn's as 
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the Pn's play regarding the pn's (see Table 1); this means that they could be obtained from a functional 
£'(.\.) = £i.\B^'.). 

Note, however, that since the matrices L and A are not uniquely defined in the d2 x d2 principal minor (in 
the Theorem we fixed the ambiguity by setting the corresponding block of L to zero), the above possibility is 
not the sole choice. 

Moreover, if some of the hard-edges zeroes of B2{y) belong to the real axis (i.e. if there are hard-edges 
onM.) then such choice is not viable because the integral defining C would be divergent at the hard-edges. In 
this case this simply means that pn 's do not belong to a biorthogonal pair for some semiclassical functional 
but are just defined (for n > (72 j by eq. (2-25) and the genericity issue cannot be resolved easily. 

Vice-versa, in the case none of the hard-edge zeroes of B2 belong to R then C indeed exists and is a 
semiclassical functional (with V2 replaced by V2 + lni?2j. The existence of the corresponding BOPs (hence 
the verification of the genericity assumption) then follows again from the result in [12]. 

3 Adjoint differential equations and the bilinear concomitant 

In this section we recall some results which -although simple- I was not able to find in the literature. We 
consider a nth order differential equations of the form 



A{8,)-xB{d^)jf{x)=0 (3-1) 

where A{D) and B{D) are polynomials and n = max(dcg(A), dcg(_B)): the reader should keep in mind the 
polynomials Ai,Bi of our matrix model. If we look for solutions written as "Fourier-Laplace" transforms 

/r(x):= J^d2/c-^-^(^) , (3-2) 

-where the contour of integration is so far unspecified-, formal manipulations involving integration by parts 
show that 

We point out that the relation between V and A,B in these formulae is exactly the same as the relations 
between the V^'s and Ai, B.i's of the first part of the paper. 
In the situation of interest to us we will have A, B reduced 

Definition 3.1 Two polynomials A, B are called reduced // the only zeroes that they share (if any) are 
amongst the simple zeroes of B. 

Lemma 3.1 Two polynomials A, B are reduced if and only if A + B' and B are. 

Proof Suppose A, B are reduced. If c is a common zero of A and B (hence simple for B) then B'{c) ^ 0: 
therefore, A{c) ±B'{c) ^ (because A{c) = and B'{c) 7^ 0). So A + B' and B do not share this particular 
zero. 
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Now let c be a common zero of A + B' and B: if it were not a simple zero of B then B'{c) = and hence 
also A{c) = 0. But this contradicts that A and B are reduced because they share a zero which is non simple 
for B. 

Viceversa: suppose A := A + B' and B are reduced. Then by the above A + B' = A and B are reduced. 
Q.E.D. 

This "duality" of the notion of reducedness will be important when considering the adjoint differential 
operator. 

We now remark that V is a rational function with poles at a subset of the zeroes of B 

r r 

B(y) =cYl{y- bjr^+' , c # , degi? = m, , m, e N . (3-4) 
j=i j=i 

d ™j , 

ny)-l]w- E T.7^A^ (3-5) 



d 



(3-6) 



W{y) := e-^^y^ (3-7) 
d := dcg(^) - dcg(B) . (3-8) 

[Here it is understood that if deg{A) < dcg{B) then the first sum in V' is absent.] 

Some of the zeroes of B{ij) may appear also as zeroes of A{y) + B'{y) and hence in the partial fraction 
expansion of V' those points do not appear. Since A, B are reduced, all multiple zeroes of B are not shared 
with A + B' . We will call the zeroes of B which are common with A + B' the hard-edge points (note that 
not all simple zeroes of B are hard-edge points, but all hard-edge points are simple zeroes). 

We now define some sectors Sj^\ j = 1, . . .pi, k = 0, . . .rrij — 1. around the multiple zeroes of B [hj for 
which rrij > 0) in such a way that 

^{V{y)) +(». (3-9) 

y ^bj, 

The number of sectors for each pole is the degree of that pole in the exponential part of W{x), that is d -I- 1 
for the pole at infinity and gj for the j-th pole. Explicitly 

d + i <^'^s(y)+ ^ + 1 < + 1 

2 < arg(y - b,) + < 2 ^ (3_jo) 







:6 C; 


k = 


0. 


.d; 


- ■ 






fc = 0, . 


. , m 



1, jeJ. 

These sectors are defined precisely in such a way that approaching any of the essential singularities (i.e. an bj 
such that rrij > 0) the function W{y) tends to zero faster than any power. 
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3.1 Definition of the contours 

The contours we are going to define are precisely the type of contours Tx.^,Ty,^ entering the definition of 
the bimoment functional C Let A,B be reduced: we then define n = niax(deg(A), deg(i?)) contours. The 
definition of the contours follows directly [3, 20]. We first remark that the weight W{y) is -in general- 
multi-valued since it contains powers like (y — cf with non-integer t; the multivaluedness is multiplicative and 
in fact is not very important which branch one chooses in the definition of the integrals (3-2) since different 
choices correspond to multiplying the same function by a nonzero constant. Nonetheless it will be convenient 
at some point to have a reference normalization for the integrals and hence we define some cuts so as to have 
a simply connected domain where W(y) is single-valued. We do so by removing semi-infinite arcs extending 
from each branch-point of W{y) to infinity: for convenience we choose the cuts approaching each singularity 
in one of the sectors, for example Sq \ and approaching infinity within S'q^ If deg(^) ^ deg(i?) — 1 then no 
sector is defined at oo and then we just choose arbitrarily an asymptotic direction for these cuts. Note that if 
deg(A) ^ dcg(i?) — 2 then the sum of the finite residues of V'Ay is zero, hence we could define the cuts as 
finite arcs joining in a chain the finite branch-points of W{y): the resulting domain is not simply connected, 
however W{y) is single valued in such domain precisely because of the vanishing of the sum of the residues 
of its logarithmic derivative. We will denote by T) the connected domain obtained after such surgery. 

In the following our primary focus is on the case dcg(A) ^ dcg(i3) -I- 1 and we leave to the reader to check 
the literature [20] for the remaining cases (only minor modifications are needed). 

1. For any zero hj of B for which there is no essential singularity in W we have two cases 

(a) If bj is a branch point (i.e. toj e C\Z) we take a loop (referred to as a lasso) starting at infinity 
in some fixed sector (e.g. Sq^'') encircling the singularity and going back to infinity in the same 
sector. 

(b) If bj is a pole of W (i.e. tg.j £ { — 1, —2, —3, . . .}) then we take a small circle around it. 

(c) If bj is a regular point ( namely toj 6 {0,1,2,...}) we take a line joining bj to infinity and 
approaching co in the same sector S'q^' as before (this case includes the hard-edge points for which 
we may say that toj = 0). 

2. For any multiple zero bj for which there is an essential singularity (i.e. for which nij > 0) we define 
rrij contours (which we call the petals) starting from bj in the sector Sq^ and returning to bj in the 
next (counterclockwise) sector. Finally we join the singularity bj to cx) by a path (called the stem) 
approaching oo within the sector S^'' chosen at point 1(a). 

3. If deg{A) ^ deg(i?) -|- 1 we define bo := oo and we take d := deg{A) — deg{B) contours starting at Xq 
in the sector 5^"^ and returning at Xq in the sector Sj^^^. 

The reasons for the "floral" names should be clear by looking at an example like the one in Fig. 1. Cauchy's 
theorem grants us large freedom in the choices of such contours; we use this freedom so that the contours do 
not intersect each other in C\{6j}j^i ... dcg(B) and do not cross the chosen cuts. 
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Figure 1: An example of contours T and T for a pair of reduced adjoint differential operators. The 
thick contours are the admissible ones for L while the thick dashed ones are the admissible ones for 
L*. Also shown in the picture are the cuts for W{y) and W{s) (line-dotted thin lines). 
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We will refer to these contours collectively as admissible contours for the differential W{y)dy. Note 
that we have defined exactly n = max(deg(A), deg(_B)) contours. 
It is a straightforward check to see that 

frix) := j^dyc^y-^'^y^ = J^c-W(y)dy , (3-11) 

all satisfy the differential equation (3-1); in these checks one is always allowed to perform integration by parts 
discarding all boundary terms because of the properties of the contours. We leave this check to the reader. 

The content of [20] (and of the fix contained in [3]) was to show that these functions are also linearly 
independent, hence providing a basis for the solution space. 

3.2 Adjoint differential operators and the bilinear concomitant 

In general, given a n-th order linear operator with polynomial coefficients 

n 

L ■.= Y,a.j{x)di , (3-12) 

3 

its classical adjoint is defined as 

n 

L* ■.= Y,{-d.ya,{x) . (3-13) 

3 

Between the solution spaces of a pair of adjoint such operators Legendre defined a nondegenerate pairing 
called the bilinear concomitant [16]. We will show that this pairing for our class of reduced operators 
admits a natural interpretation as intersection pairing. 

We begin by noticing that in our case the pair of adjoint operators is written 

L := A{d^) - xB{d,) , L* := A{-d^) - B{-d^)x . (3-14) 

Since A,B are reduced then L* is also reduced since 

L* = A{-d^,) - - xB{-8^) (3-15) 

in view of Lemma 3.1 (here the polynomials are A{—y) — B'{—y) and B{—y) which are clearly reduced iff 
A{z) — B'(z) and B{z) are). Therefore L* is in the same class of operators as L and can be solved by contour 
integrals in the same way. The solutions of L*g = are of the form 



(3-16) 



y('y--=§^^=v'i')-i^^B(.s)y. (3-17) 

An inspection shows that the sectors around the multiple zeroes of B{s) where ^{V{s)) -co are precisely 
the complementary sectors defined in (3-10) for V. We normalize V{s) by choosing the integration constant 
in such a way that 

W(s) := e^("^ = ^e^("l (3-18) 
B{s) 
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(here is supposed to be defined on the simply connected domain P). One then proceeds in the definition 
of the admissible contours f for the weight W{s) and of the simply connected domain V in exactly the same 
way used for W{y). We make the following important remarks: 

1. If bj is a hard-edge point for W{y) (i.e. it is a zero of B{y) but a regular point for W{y) where W does 
not vanish) then bj is a simple pole of W{s). 

2. If bj is a zero of multiplicity m of W (y) (i.e. a simple zero of B{y) such that the residue of {A + B')/Bdy 
is a negative integer) then it is a pole of order to + 1 for W{s). 

3. In all other cases, the type of singularity of W and W is the same (logarithmic branch-points or essential 
singularities of the same exponential type). 

4. The intersection T) n V \s the disjoint union of simply connected domains where W{y)W{y)B(y) is 
constant. These constants depend only on the residues of V'{y)dy mod Z. 

These observations and the fact that B {y)W {y)W {y) is locally constant (where they are both defined) follows 
immediately from their definition and eq. (3-18). 

From the definitions of the contours it is not difficult to realize that dual contours can be chosen such that 

1. For each flower (petal + stem) one can choose a dual flower whose elements intersect only the arcs of 
the given flower. (This includes the petals at oo, in the case dcg(^) ^ dcg(i3) -I- 1). 

2. For each pole c of W{y) (whose corresponding admissible contours F is a small circle) the dual admissible 
contour for W{s) is a semi-infinite arc starting at c and going to oo and can be chosen so that it intersects 
only its dual. 

3. For each zero or hard-edge point a of W{y) (whose corresponding admissible contour is a semi-infinite 
arc starting at a) the dual admissible contours for W{s) (which is a small circle around a) intersects 
only r. 

4. For each non-essential other singularity of W{y) (i.e. a simple zero c of B{y) such that the residue of 
{A + B')/Bdy is in C\Z), where the admissible contour F is a lasso around c, the dual loop F (also a 
lasso around c) is also chosen so that it intersects only the dual lasso (at two points). 

Lemma 3.2 Consider the two adjoint differential equations 

'A{d^)-xB{d:,)Jf{x)=0 (3-19) 
^(-5,) - B{-d^)x] g{x) = . (3-20) 
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The solutions are of the forrri^ 



-,-V{y)+xy 



dy, V 



J 



A{y)+B'{y) 



dy 



(3-21) 
(3-22) 



;r J B{s) 

Then the following expression is constant and defines a nondegenerate bilinear pairing (the bilinear con- 
comitant ) between the solutions spaces of the two adjoint equations: 



B(y) - Bis) 



y-s {y -sf 



A{y) - A{s) - B\s) 



y — s 



x(y-s)-V(y)+V{s 



Proof. The integral representation of the solution is easily verified. We now write 

= g{x) ixB{y) - A{y)) e-^('')+^M2/ 
= f{x) j {xB{s) - A{s) - B'{s)) e^(")-^Ms 
We take the difference and obtain 

= J (x{B{y) - B{s)) - iA{y) - B'{s) - A{s))) e^^^-^'^-^^^'+^^^Myda 



)d2/d^-23) 

(3-24) 
(3-25) 

(3-26) 



It is promptly seen that the integrand of this double integral is absolutely summable w.r.t. the arc-length 
parameters along F and F, hence we can integrate w.r.t. x under the integral sign, thus obtaining the bilinear 
concomitant; 



{B(y)-B{s)) 



1 



y-s {y-s)' 



A{y)- B'{s)- A{s) \ ^x(y-s)-v(v)+v(s] 



y-s 



'dsdy (3-27) 



Note that the expression under integration is regular at y = s, and is -in fact- a polynomial in ?/, s 



{B{y)-B{s)) 



1 



y 



{y - s)^ 



A(y) - B'{s) - A(s) 
y-s 



- xB'is) - Ib"{s) - A'{s) + 0{y - s) 



In particular the integrand is absolutely integrable w.r.t. the arc-length parameters and hence the order of 
integrations is irrelevant. This concludes the proof. Q.E.D. 

The bilinear concomitant is -in a certain sense- an integral representation of the intersection pairing of 
the contours of integration. To make this statement more precise we first prove the following standard 

Lemma 3.3 Let Q{y,s) be a meromorphic function V x V where V and V are simply connected domains 
and with the only singularities being a double pole as y ^ s (in T> "D). Suppose that in each connected 
component ofVr^T) there is a constant c such that 

c 



niy,s) 



(y-s) 



+ 0(1) 



(3-28) 



'^The formula depends on the integration constant in V, namely these solutions are defined up to multiplicative constants 
since they are solutions of a homogeneous linear ODE. 
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as y ^ s within ttie intersection domain. Let T cz V be a smootli curve such that 

n{y, s)dy = (3-29) 



I 



Jr 

Let T <z V be a curve of finite length intersecting once T at p and oriented positively w.r.t. T: then 

J dy ^ dsf7(y, s) = 2inc{p) (3-30) 

Proof. The integral 

/(s):=J^fidy (3-31) 

defines -in principle- different holomorphic functions in the connected components of P\r: the difference 
among them -however- is the residue 

res n{y, s)dy (3-32) 

which is zero by the assumption on fl. Hence the analytic continuations of /(s) from one component to the 
other all coincide. In our case they are all zero. The key fact is that, since fl is singular on the diagonal, the 
orders of integration matters (otherwise (3-30) would give zero by interchanging the order of integration). 
We compute the integral as a limit of regular integrals where we can interchange the order of integration 

(3-30) = lim f dy r dsf7(y,s) , (3-33) 

where is the curve (or union of curves) obtained by removing a small e-arc (which we denote by F^, i.e. 
an arc from p — e to p + e, where these two points lie on the curve F at distance |e| from the intersection 
and the direction of e is the same as the orientation of F) around the intersection point p. This allows us to 
interchange the order of integration under the limit sign 

lim dy dsfi(y, s) = lim ds dyfl{y,s) = — lim ds dyQ(y,s) = 

where we have dropped the 0(1) part since the length of F is finite and that of F^ tends to zero. In the last 
expression the inner integral is -strictly speaking- defined only for s # p: however on the "left" and "right" 
the result is the same and gives 



r r c(p) . f 

lim ds dy- -:r = c(v) lim ds 



b - p + e 



c{p) lim In ('^^^1 - In ( '-^1±1 ] (3-35) 



\a — p — e J \a — p + e 

In this last limit the logarithms appearing have different branches: in particular the second differ by 2i7r from 
the first, hence the result follows by taking the limit. Q.E.D. 
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We now come back to the computation of the concomitant: first of all, since we know that the result is 
independent of x we set a; = 0, so that we have to compute 



B{y)-B{s) A{y)-A{s)-B'{s) 



(y - 



y 



^dyds 



(3-36) 



We have already remarked that this integral can be computed in either orders and gives the same result. We 
express it in terms of 



ds 



Bif,9) = (2)-(1) 

(2) := J^dyJ^ds 



B(y) A(y) 



(y - s) 

Bis) 



y-s 

A{s)+B'{s) 
y-s 



W{y)W{s) 

W{y)W{s) 



(3-37) 
(3-38) 

(3-39) 



The integral (2) is zero because the inner integral w.r.t. s defines (for y ^V) the identically zero function, as 
it is easily seen after an integration by parts. The integral (1) is computed using Lemma 3.3 after noticing 
that 



B{;y) A{y) 



{y-sy- 



y 



W{y)W{s) 



B{s)W{s)W{s) 



+ 0(1) 



(3-40) 



{y - s)2 

and hence satisfies the condition of the Lemma for fl. The contour F satisfies the condition of the Lemma. 
The contour F is not necessarily of finite length, but we can take only a small arc around the point of inter- 
section and the remainder will be computed to be zero by interchanging the order of the integrals. To rigor 
one should also consider the common endpoints of contours like the petals and dual petals: it is easily seen, 
however that those points do not correspond to a singularity of the integrals (w.r.t. the arclength parameters) 
because of the fast decay of the weights W and W. For example, if the two contours F,F form an angle 
9 e [0 -I- e, TT — e] (asymptotically) near a point b (where W, W have an essential singularity) then 



W{y)W{s) 



iy 



W{y)W{s) 
^ — [see fig. 2]. 



sin"' 6*1?; - fej^ 



(3-41) 



which is still jointly integrable w.r.t. the arc lengths (recall that the directions of approach of F and F are 
such that the weights tend to zero faster than any power of the local coordinate). 

It is then clear that if F F are a circle and a semi-infinite arc (or vice-versa) the bilinear concomitant for 
the corresponding dual solutions is a nonzero constant (which depends on the choices of the branches of W 
and W). This is immediate for a pair of contours which intersect only once. For a pair of lassoes (which 
intersect twice and with opposite orientations), calling pi,p2 the points of intersection we have 



^(/r,.9f) = ±(^^(^1)1^(^1)5(^1) - W{p2)W{p2)B{p2)) 



(3-42) 



Since the local behavior at the singularity embraced by the lassoes is a noninteger power, let's say {y — cf , 
then the values of BWW on the two intersection points (which lie on different sizes of the union of the cuts 
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s 



Figure 2: Illustration for estimate (3-41) 



for and W) satisfies 



Wipi)Wipi)Bipi) = e^"'W{p2)W{p2)B{p2) 



(3-43) 



so that 



S(/r,5p) = ±(M/(pi)Ty(pi)i3(pi)(l -c^"*) ^ 



(3-44) 



For dual flowers it is convenient to choose different paths for the dual contours as shown in Fig. 3, where the 
petals have been replaced by stems using a linear combination of the contour-integrals of the same petals and 
stem. The sub-block of the concomitant involving these contours is nondegenerate, since it can be given a 
diagonal form with nonzero entries on the diagonal. The precise values are not important since we are free to 
re-scale each solution /r and gr- Summarizing we have proved that 

Proposition 3.1 There is a normalization of tlie integrals /r and such that the bilinear concomitant is 
precisely the intersection pairing of the contours T andT. With appropriate choice and labeling of the contours 
the pairing is represented by the identity matrix. 

Remark 3.1 The content of Prop. 3.1 is that if the solutions /r and correspond to contours that can 
be deformed (by Cauchy's theorem and without changing the analyticity properties of the functions /r,5f 
respectively) In such a way that they do not intersect, then the bilinear concomitant of this pair is zero. 

Viceversa, if this cannot be done, the bilinear concomitant is nonzero; we can always choose the contours 
and dual contours in such a way that each contour intersects one and only one dual contour. For example the 
equivalent choice of contours to Fig. 1 is given by the arrangement in Fig. 4. 

Note that two dual lassoes intersect at two points but that -by virtue of (3-44) their "weighted" intersection 
number is nonzero (whereas the usually defined intersection number would be zero). 

4 Auxiliary wave vectors 

Caveat In this section we will make statements concerning the biorthogonal polynomials p„, s„ and the 
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Figure 4: The arrangement of dual contours for the same example as in Fig. 1: in evidenee only the 
different choice of admissible dual contours at the "flower" . 
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corresponding quasipolynomials i/)„,0„. It will be understood that 

1. Any statement made on the t/v/s and the Fourier-Laplace transforms of the 0„'s admits a specular 
statement for the (p^'s and the F-L transforms of the V'n's- 

2. Any statement made on the i/'n's admits an analog statement for the V'n's and ipn's by replacing the 
moment functional C with £ or £, and specular statements for (j)n,4>n- 

Consider the functions 

BA^;y,s) :^ f^^l^^ ^ V Mvl^^ifl^AM) (4-1) 

V y-s \ y-sj y-s ) 

i^P ■■= ^ I ds JJ d^dyB^ix; y, ,),^y-^^s-v.(yH%(s)^ (4_2) 

If X belongs to a contour F^ ^ of the integration we obtain 

^l!\x)+ = +;^62(f,r^,,)x^,,7/;„(x) (4-3) 

where the subscript x+ denotes the boundary values from the left/ right and S2(r, i,) stands for the constant 
(in a;) bilinear concomitant 

S2(f,r^,,) := ^ [ds r dyB2(x;y,s)e^^(«)-^^fe)+-fe-«) . (4-4) 

Therefore their jump across the contours of discontinuity is a constant multiple of ipn{x). 
We have 

Proposition 4.1 The sequences of functions {ipP^}neN satisfy the same recurrence relations (for n large 
enough ) as the quasipolynomials ipn 

(12 ^ \ d2 

+ 2 (n)^i^_^ J = J] a, (n)^fj^. , n ^ + q2 (4-5) 

8. (^i^P + £ ™, (n + j>f+^ j =-Y^p^(n + j>f+^ , n^l (4-6) 
(For the proof see App. A.l) 

Definition 4.1 Beside the wave vector we define the following d2 auxiliary wave-vectors 

^^"^2^) := 7^ f ds rr dedyB2(x;2/,s)c«^-"*-^=^(«)+^^W^-^(0 , u = l,...,d2 (4-7) 

^(°)(2;) := ^(x) . (4-8) 
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We also define the dual wave vectors 



:=e^i(^' ff c^^-^^'^^^ ^—<i.(y)d^dy (4-9) 
t(^)(x):=r dye'^y^iy) ,iy = l,...,d2 (4-10) 



Proposition 4.2 T/je components of the dual wave vectors satisfy the recurrence relations 



X 



l''^ + f]^An+j)ti'l)= I; S,(n+j>M^.+,5.o'5„oV^e^^(^) , = 0, . . . , (4-11) 

+ £ "^^w^i-, J = £ /^.w^":!, = i,...,d2 . (4-12) 

Remark 4.1 The wave vector c/oes not satisfy a finite-term differential recurrence relation: a formula 

X 

can be derived but it is not useful for our purposes. 

Proof The formulas for the Fourier-Laplace transforms follow from integration by parts from the relations 
satisfied by (pniy) (Prop. 2.3). We only point out that integration by parts does not give any boundary 
contribution because s„ -I- '^ij{n + j)sn+j{y) is divisible by B2{y) and hence vanishes at the hard-edge 
end-points. 

The only relation that needs to be checked is the multiplicative relation for v = Q. Denoting temporarily 
by a tilde the linear combination 



bn+Yj^3{n+ j)(j)n+j , (4-13) 



we have 



=e^^(")JJ e«^-^^(«)^„(y)d^dy + e^i(^) JJ e«''-^i(«)-^^„(2/)ded2/ = 

=e^^<")5«oV^+ 2 + ■ (4-14) 

4.1 Christoffel— Darboux identities 

In the general theory of the two-matrix model the following kernel plays an essential role in the computation 
of statistical correlation functions 

K^2{^.y) ■■= 2 P,(a;)5,(y)e-^^(-)-^^(^) = ^ ^A^mv) ■ (4-15) 
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In a previous paper by the author and collaborators [4, 8] the case of polynomial potentials Vi was considered 
(without hard-edges) and it was of capital importance the existence of a Christoffel-Darboux identity allowing 
to express /vf^ (or rather some transform of it) in terms of bilinear combinations of the BOPs involving only 
a number of BOPs depending only on the degrees of the potentials. 
We look for a similar bilinear expression in this model. 

Definition 4.2 We define the windows of the wave vectors v]/(^) and ^ = 0, . . . ^ d2 

"z a 

$(r'(x) := [^l^l, . . , ^i^Hx) := [V^^),^, . . . ,^('')]* . (4-16) 

We rewrite (4-15) in terms of the wave vectors 

Ar, = $'(y)n^.vl.(x), U^:=\ ' (4_i7) 
X cc I otherwise . 

Recall the multiplicative and differential recurrence relations in Prop. 2.1 and Prop. 2.3 (which we rewrite 
here for the reader's convenience) 

5y$*(l -f£) = , x{l + L)\if = Aqi 

cXj 'X X X 

(l+L)-iA = A(l+£)-i =:Q . 

Consider now the following expressions 



{x + dy)<^\y){l + L)n(l + L)-i*(x) 



y 

X 



=$*(?/)(l + L)n(l + L)-M(l + L)-^^{x) - $*(y)^n(l + L)-^^{x) = 

'X -Ji X X 

=$*ni(i + L)-i* + $[£, n](i £)-^A(i l)-i*-f 

X X, X 'X 

-$*nl(i £)"^^' - n](i £)"^^' = 

X X X X 

=$[£,n]Q(i + L)-i^ - n](i + L)-i^ = 

X X X X 

=$[£,n]Q$-$*[A,n]$ (4-18) 



where we have set Q := (1 -|- L) ^A. We now use the fact that Q is the recurrence matrix for the associated 
vj/ wave vector (see Prop. 2.2 where vj/ := pe~^^(^^) and obtain 

X' X 

{x + dy)<^'{y){l + L)U{1 + L)-i*(2;) = <^{y)[xL - A, U]^{x) (4-19) 



$ = (1 L)-^xiJ = p(a;)e-^i(^) (4-20) 

X X, 



X X 

With these notation we have 



$ = (1 L*)$ = s{y)B2e-^^^y' = s{y)e-^^^y> (4-21) 

X X 
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Theorem 4.1 (Christoffel-Darboux identity) For the kernels 

^'"{x)Un^{x') , (4-22) 
{xfuN^ix') , J = 1, . . . , d2 (4-23) 

= $(^^(a;')*A(x)§(a;) (4-24) 
= $(^'(x')*A(x')$(a;) . (4-25) 

fnofe the argument of An in the two formulae) where An{x) := — xL, IIn~^ . 

Proof The identity for K^{-' {x,x') follows by performing integration by parts on (4-19) and noticing that 
the boundary contributions vanish since ^{y) = i?2(y)s(y)e^^^^^' and B2{y) vanishes at the hard-edges. The 
identity for K^^-' (x, x') follows from the one for K^{'' and this manipulation 

(x' - x)^^'Hx'Y{i + L)n§ix) = 

{x' -x)(^^^\x'Y[L,Il]^i{x) +^^^\x'yum{x)] = 

\x ^' X J 

{x' - x)K^{'{x',x) + {x' - x)^^'\x'y[L,U]^i{x) (4-26) 

00 CO 



(x'-x)$(^''*(a;')[£,n]$(x) =$(^)'*(a;')A(a;')$(x) (4-27) 
Q.E.D. 

Note that -with a slight abuse of notation- in the RHS of the CDIs we can replace the wave vectors $ by 

'X 

the corresponding window since the matrix A,i has a nonzero square block of size d2 + 1 with top-right 
corner in the (n — l,n) entry, and hence the bilinear expression $A\I/ only involves the terms in the dual 

windows and 'i'n- We will denote from now on by A only the d2 + 1 square matrix which is relevant to 
the pairing. 

The importance of the theorem is that we can express the kernel Ku in terms of the dual quantities (f'j^i^) 
and ipnix') involving only the indexes N — d2 ^ n ^ N. 

Note, however, that we must introduce the orthogonal polynomials p for the associated moment functional 
C in order to find a Christoffel-Darboux relation similar to the standard one for orthogonal polynomials. 



K^{''{x,x') = e-y$\y)UN^{x') = 
K^^''{x,x')= f e^y^'{y)UN^{x') = 
we have the identities 

{x' -x)K^{^x',x) = 
{x' -x)K^{\x',x) = 



{x' -x)%^'\x'YU^{x) 



so that 



{x' -x)K^{^{x',x) = 
= {x' -x)K^{'(x',x) 
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Theorem 4.2 (Auxiliary GDIs) The auxiliary wave vectors enter in the following auxiliary ChristofFel-Darboux 
identities 



(a) (z-x)$(°)(z)*n„*(o)(x) =$i°)(z)A(z)§„(a;) +e^i(")-^i(") 
(z -x)i^°\z)*n„$^°'(x) = $i°)(z)A(x)§„(x) +e^i(^)-^i(") 

(X, « 

(b) (z-x)$«(z)*n„M.('=)(x) =ti^'(z)A(z)§„(x)--^ r r S2(x;y, 

(z -a;)l''\z)*n„$^''(a;) = $i^)(z)A(x)$„(x) - f f 62(2:; y, 

j,k = l,...,d2 . 



(4-28) 



yz-a;s + V2(s)-V2(y) 



s e 



(4-29) 



(For the proof see App. A. 2). 
4.2 Ladder matrices 

In this section we derive an expression for the ODE satisfied by the polynomials in terms of the so-called 
"folding" (see [4]). This will have certain advantages when explaining the relations between the various ODEs 
that naturally appear in the problem: a different explicit representation of the ODE will be given in the next 
section as well, using a completely different approach based upon the explicit integral representations of the 
wave vectors and on the duality provided by the Christoffel-Darboux pairing. 
We first have the simple lemma 

Lemma 4.1 (Ladder matrices) The multiplicative recurrence relations for the wave vectors xl/C^),^ = 

<X rx rxj 

(J=l,...d2) 



V{1 + L)^^J(°'^ = yl^f'") 



:(1 -fl*)$(^) = 



are equivalent to the relations 



^<fH^)=tif|i(a:)aJx) 



(4-30) 



(4-31) 
(4-32) 



where 



a„(a;) = A ■ 



1 

q;_i(?i) 



[ad2 



,ao{n)] + 



o:-i{n) 
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in (a;) = A- 



a_i(n — 1) 



So(n) 
ai(n-l-l) 

a<i2 {« + 'i2) 



[1,0,...,0] + 



a_i(n — 1) 



^i(n + l) 




[l,0,...(i]34) 



and A denotes the upper shift matrix (of size d2 + l). The relations (4-31 ) and (4-32) hold also for the other 
sequences of windows ^E'l''' and provided that n ^ d2 -\- q2 (n ^ I respectively). 

Proof. The proof follows immediately from the recurrence relations for the wave vectors vj/'"' (the quasipoly- 

X 

nomials) and ^"^^^ (the Fourier-Laplace transforms) by solving for ■(/;„+i(x) (or0^^ ^) in terms of ?/'n-d2: ■ • ■ jV'n 
{(j)^, . . . , 0^^^^ ) and rewriting the relation in matrix form. The statement for the other sequences of windows 
follows from the fact that the corresponding wave vectors satisfy the same finite-term recurrence relations in 
the specified range (see Prop. 4.1 and Prop. 4.2). Q.E.D. 

Lemma 4.2 (Folded recursion relations) The differential recurrence relations for the wave-vector ^ 



are equivalent to the relations 



Mn ■.= 1 + Yj mj{n)sLn ■ ■■Sin+j-l 

rn.j(n) := dia.g{ifij{n -\- j - ^2), . . . ,m-i{n -\- j)) 
Bn ■■= ^_i(n)(a„_i)"^ + \{n) + ^ ^j(n)a„ • • •a„+j_i 
:= diag(/?j (n j - ^2), ■ • .Jj{n + j)) 



(4-35) 

(4-36) 

(4-37) 
(4-38) 
(4-39) 

(4-40) 



Proof. The formula is an iterated application of the ladder recurrence relations (on a window of consecutive 
elements with indexes n — d2, - ■ ■ ,n) to the differential recurrence relation for the wave vector (see [4] for 
more details). Q.E.D. 

Remark 4.2 A completely analogous statement can be derived for the windows of the dual vector 
j = l,...,d2. 

Remark 4.3 The matrices a„ have a companion-form and are invertible since the determinant is — 2''\^(n) 
which has been proved nonvanishing in Thm. 2.1. Moreover the inverse is also linear in x (the details are left 
to the reader). 
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Remark 4.4 By the very definition Mnix)'^r, 
functions) for tfie moment functional C. 



is the window of quasipolynomials (and associated 



Corollary 4.1 The d2 + 1 columns provided by the windows of the auxiliary wave vectors vl/(j)(a;) provide a 

CO 

fundamental system for the ODE (4-36) for n ^ d2 + q2- 

Proof. From Prop. 4.1 we know that the components of the auxiliary wave vectors satisfy the same recurrence 
relations (both multiplicative and differential) as the quasipolynomials provided n is large enough. Moreover 
the recurrence relations always involve a fixed number of terms with indexes "around n" : since the derivation 
of the ODE is entirely based on the recurrence relations the statement follows. Q.E.D. 

Proposition 4.3 The determinant of A4n{x) is proportional to Bi{x) by a nonzero constant. 

Proof. Consider the window of polynomials p„ := [pn-d^i ■ ■ ■ iPnY'- from the definition of the matrix M it 
follows that 

Mn{x)l>n{x) = Bi{x)Yin{x) (4-41) 

We first prove that dct7Vl„ (which is a fortiori a polynomial) is divisible by Bi. Let c be a zero of Bi 
of multiplicity r: at least one component (say the ^-th) of Pn(c) is nonzero because of the very genericity 
assumption which guarantees the existence of M (2-38). Let E{x) be the matrix obtained by replacing the 
£-th column of the identity with p„(x). Clearly dctE(x) is nonzero in a neighborhood of a; = c by our 
definition of £. It follows that the £-th column of M.nE is precisely SiPn and hence each component vanishes 
at c of order r. Also 

det MnE = pn-d2+e-i{x) det A4„ (4-42) 

and pn-d2-i+e{c) 0. On the other hand dct^AnE must vanish at x = c of order r since the whole i-th 
column does. Repeating this for all roots of Bi we find the assertion of divisibility of det A^„ by Bi{x). 
On the other hand, using a technique of evaluation of determinants used in [4], 



det Mn = det 



L(d2-l-l)(9l-l-l) 



























rrig, (n) 




mi{n) 






(4-43) 



Considering carefully the structure of the sparse matrix in the last identity, one realizes that the highest power 
in X is 

detAl„ = + + 0(^91-1) (4_44) 

This shows that (since the coefficient does not vanish as per (2-19,2-48)) the determinant is of degree 
qi = degi?i; since it must be also divisible by Bi, this concludes the proof. Q.E.D. 
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Corollary 4.2 The windows 'i'n, satisfy 

5,*„ = -yW-i (^Bn + dMn^ (4-45) 

d^^n = -BnM-'^n (4-46) 

where Bn,Mn are defined in ((4-37)-(4-40)). The ODEs have the same singularity structure as V{. 

The first relation follows from (4-36) and the second from the fact that A^„(a;)^„(a:;) = ^'„(a;). 

This shows that the ODE's for 'i'n and are gauge-equivalent, the gauge being provided by the (poly- 
nomial) matrix A4n- Moreover formula (4-46) together with Prop. 4.3 shows that the singularities of the 
differential equation are at the zeroes of Bi{x). 

4.3 Differential equations for the dual pair of systems 

In this section we present an explicit formula for the ODE satisfied by the dual pair of fundamental systems, 
in particular the polynomials ^/'„ and the Fourier-Laplace transforms 4>^'^- The result generalizes those of [7] 
but the method of derivation is similar to the one adopted in [8], with additional complications deriving from 
the presence of boundary contributions in the integration by parts at various steps of the derivation. 

Notation. In the proof of this and the following theorems we will encounter integrations by parts that yield 
nonzero boundary contributions. Typically we will encounter integrals of the form 

JJ ye'"'-^^('')F(p)0„(y)d2/dp , (4-47) 

where F{p) is some expression (typically polynomial or rational in p) possibly depending on "external" variables. 
If we attempt an integration by parts on the term ye'^^ = dpe^^, we obtain a certain number of boundary 
terms. In all cases they will be boundary evaluation on the various contours r^.^; it is the nature of all these 
integrals that only the contours emanating from a hard— edge point give a contribution, due to the fast decay 
of c^^'i*'') at all the boundary points of the other contours. In the above example and in all minute detail, we 
have 

JJ yefy-'''^P^F{p)(j>^{y) = - JJ ePy-'^'^P^ {-dp + Vl{p))F{p)<j>,n{y) + (Boundary terms) 



di 6.2 o 

(Boundary terms) = ^ c-^i('')F(p) x^^^ I 
/i=i v=i ''r 



(4-48) 



The evaluation at the boundary points of the various contours yx,^l is clearly to be understood as limits 
along the contours; the decay of c"^^^'') along the contours gives zero contributions except for the hard-edge 
contours, at the (finite) boundary of which Vi{p) is regular. In order to economize on space, we introduce the 
following shorthand notation for the above boundary terms 



:= (Boundary terms) (4-49) 
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Theorem 4.3 The dual fundamental system. 



r 


.. 




Vn-l 




■ Vn+d2-l 








Vn-1 




■ yn+d2-l 






Ad2) 


L 







(4-50) 



satisfies the ODE 



V(ix) 

Pn.n — 1 Pn.^i 







■ ■ ■ -fn,n + c/2 — 1 

Pn + l,n '■ 



Pn + d, 



+ diag{Pn + d2.n-l, ■ ■ ■ , Pn + d2,n+d2-l)^n ^ i^) + A(x) 



X — p 



Wab{x) ■= C[ p„-d2+a{p) 



p — X 



p€<:\c>r: 

s„_i+h(y) 1 , a, 6 = 0, 1, . . . ,d2 



Pj.k ■■= ^Pj\ysk) ■ 

where is the ladder matrix for the dual wave vector (Note that P = {{1 + M)^^BY) 
(For the proof see App. A. 3). 
Theorem 4.4 The direct system 



A{x)W(x) 

(4-51) 
(4-52) 



satisfies the ODE 



Vn-d2 



Vn-d2 



2(d2) 
Vn-d2 



(4-53) 



Pn — d2.n—d2 
-P7i — d2 + l-n—d2 




Pj,fc := C{pj\ysk) , 



• • - Pn — d2:n — l 

: 

Pn-l,n-2 Pn-l,n-l 

PrK«-l V{{X) 

A{x) + W(x)A{x) 



+ diag(P„+i 

,n-d2 1 • ■ • I Pn + l,n)^n-l~^ 



(4-54) 
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where W{x) was defined in ttie previous ttieorem and a„_i is the ladder matrix implementing the multiplicative 
recurrence relations = a„_i5',i_i as per Lemma 4.1 (in particular eq. (4-31)) specified to the hat-wave 
vectors. 

(For the proof see App. A. 4). 

5 Dual Riemann— Hilbert problems 

The shape of the Christoffel-Darboux identity (Thm. 4.1) suggests that the duality of the Riemann- 
Hilbert problems (and of the differential equations) involves naturally the dual pair of fundamental systems 
^„(a;), 'S'„(.t) defined in Thm. 4.3 and Thm. 4.4. Recall (from Section 3) that we can choose a basis in the 
relative homology of contours Ty^^ and (and a rescaling of the $ wave vectors depending only on the 

'X' 

residues of ^2 (y)dy) which span the solution space of the two adjoint equations and with bilinear concomitant 

l32(Xy,^,^y,p.) ■■= ^y,4'^y,^L = 5 ^iy ■ (5-1) 

We can rewrite (Thm. 4.1) as 

n-l 

{x - x') 2 0^ (^') = ^ (^)A(x')*l°' {x') (5-2) 

i=o 

{x - x') 0W {x)i^f {x') = ^l:^ (x)A(x)$(o) {x') (5-3) 

3=0 

V = 1, . . . ,c?2, where we stress the fact that on the LHS we have the quasipolynomials ■0„ whereas on the 
RHS we have the i/'n's. 

Theorem 5.1 The fundamental dual pair is put in perfect duality by the Christoffel-Darboux matrix A 

*„(x)A„(a;)*„(x) = 

where B2{», •) represents the (constant in x) bilinear concomitant for the solutions of the adjoint ODEs along 
the contours Ty_^, Ty^^, = 1, . . . ,d2. By suitable choice of the homology classes we have seen that we 
can always assume it to be diagonal. The entries on the diagonal are nonzero and may be set to 1 by suitable 
rescaling of the d2 left-most columns of'ifn: these re-scalings depend on the way we have performed the cuts 
in the definitions of V2 and V2 but depend only on the residues ofV^ mod Z. 

(For the proof see App. A. 5). 

5.1 Riemann— Hilbert data 

In this section we summarily indicate how to obtain the data of the Riemann-Hilbert problems solved by the 
dual fundamental systems. The details are considerably involved and not strictly necessary in this paper. They 
will appear in a different publication. 



B2{;' 



(5-4) 
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Since the two matrices and are put in perfect duality by the Christoffel-Darboux pairing, it is -in 
principle- sufficient to describe the Riemann-Hilbert data of one of the two members of the pair, the data for 
its partner being completely determined by duality. 

It is significantly simpler to analyze the RH data for the matrix We recall that this means controlling 
the jump discontinuities and the asymptotic behaviors near the singularities. 

Jump discontinuities. They are uniquely due to the first row in the definition of ^„ and occur at the 
contours F^,!/: 

1 2iTr>c^^i 2i'iTK,j,2 ■ ■ ■ '^'ii^Xu.d^ 
1 



(5-5) 



where x+ denote the boundary values on the left/right of the point x e Tx,y- 

Note that the fundamental matrix ^n(x) satisfies a similar jump condition which can be read off eq. (4-3) 
(specified to the t/)„ quasipolynomials). 

Singularities The bottom d-i rows (the Fourier-Laplace transforms) are entire functions. The only 
singularities in the finite part of the plane arise from the first row apart from the jump discontinuities 

(discussed above) we have all the singularities of c^^^^) and the logarithmic branching singularities around the 
hard-edge endpoints. Note that the (piecewise analytic) function 



i 



<i>n(y)- 



-^^(-)$„(a:) 



(5-6) 



has a well defined limit as x approaches any of the non hard-edge endpoints (where it is understood that the 
approach occurs within one connected component of its domain of analyticity). Indeed, if c is such a point 
one finds 

i^„(c) = JJ €>„(y) ^ ^ (5-7) 



which is a well-defined value. In other words, near a non hard-edge singularity one has 

*„(x) ^ diag (e^i— 1, . . . , l) yo(l + 0{x - c)). 



(5-8) 



where Yq is just the evaluation of the Fourier-Laplace rows and the Fn{x) defined above at the point c, and 
Visaing denotes the singular part of Vi at c. 

Near a hard-edge point a; = a, if Tx.va is the the hard-edge contour originating from a, we find that the 
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matrix 



Y{x) 



1 ln(a; — a)K^^^ 



1 



(5-9) 



has a removable singularity at a; = a and from this we can obtain the asymptotic behavior near the hard-edge 
endpoints. 

Stokes Phenomenon. Possibly the most intricate part is the description of the Stokes' phenomenon at 

a; = 00. 

Indeed, apart from the aforementioned jump-discontinuities of in a neighborhood of cx) (which may be 
interpreted as part of the Stokes data), the first row displays no Stokes' phenomenon, and has an asymptotic 
behavior which encodes the orthogonality 



-Vi(5)-^Cy 



/^e^i(^)x-"-i(l -FO(1/.t)) 



(5-10) 



The remaining part of the Stokes phenomenon is given by the asymptotic behavior of the d2 Fourier-Laplace 
transforms: this is precisely the same Stokes' phenomenon displayed by the solutions of the ODE 



{A2(d^)-xB2{d^))f = Q 



(5-11) 



These solutions are described by contour integrals of the same kind as the ones appearing in the expressions 
for a standard steepest descent formal argument shows that the leading asymptotic is determined by the 
saddle-point equation 

A2{y)+B'2{y) 



Vi{y) = X 



(5-12) 



B2iy) 

{x —>■ oo) which has d2—H solutions {H being the number of hard-edge contours, i.e. the number of (simple) 
zeroes of i32 which cancel against corresponding zeroes of the numerator in (5-12)). 

Whereas it is not very difficult to analyze the formal properties of the asymptotic, it is considerably harder 
and outside of the intents of the present paper to present the Stokes matrices associated to this Stokes' 
phenomenon. We leave this topic to a different publication. 



5.1.1 Isomonodromic deformations 

The (generalized) 2-Toda equations for this reduction as explained in the introduction, determine the evolution 
of the biorthogonal polynomials under infinitesimal deformations of the parameters entering the semiclassical 
data Ai,Bi. It is more convenient to parametrize the polynomials Ai, not by their coefficients but by the 
location of the zeroes of Bi and the coefficients in the partial fraction expansions of the derivative potentials 
V/. Following the strategy in our [4, 9, 11] one could easily write the pertinent 2-Toda flows corresponding 
to these infinitesimal deformations. 
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At the level of the pair of fundamental systems the flows will generate isomonodromic deformations for the 
ODEs satisfied by ^„ and provided that the exponents of formal monodromy at the singularities remain 
unchanged. In this case these are precisely the residues of V({x)dx and V2{y)dy at the various singularities. 

The reason why the deformations are isomonodromic is that -by their very definition- the fundamental 
systems are functions of these deformation parameters and the matrices .^„^„~^ (and ^„^„ , the dot 
representing a derivative w.r.t. one of the monodromy-preserving parameters) are rational (or polynomial) 
functions of x, which follows from the analysis of their behavior at the various singularities ([17, 10] for details 
on the general properties of isomonodromic deformations). 

The details of this isomonodromic system could be derived from the complete Riemann-Hilbert charac- 
terization of the fundamental systems and are beyond the scope of this paper, although their derivation is -in 
principle- a straightforward computation. 

A Proofs 

in this appendix we report all proof of more technical nature. The expressions are rather long and hence 
to shorten them we have decided to suppress explicit reference to the variables of integration in the multiple 
integrals below, since which variables are integrated on which contour is unambiguously implied by the context. 
We have adhered to the following general naming scheme: the variables ^, p are integrated along the contours 
Fa; 1, appearing in the integral the variables y and i] are variables integrated on the Fy ^'s. The variable s is 
always running along the dual contours f ^,^ (the admissible contours for the differential W{s)ds = c^^'^Ms = 



A.l Proof of Proposition 4.1 

We temporarily denote by a tilde the following linear combination 

^ 92 

= Ipn +^ij{n)lpn-j (1-1) 
1 

and notice that 

xi>n = 2aj(n)7/i„_j . (1-2) 

-1 

For the transformed functions V'n"' (denoting by a tilde the same linear combination) 

^^(H = —( {{ S2(x;2/,s)c«^— = (i_3) 

r rr B,{x; y, ,)e««— ^^(^)+^^(^) ( MO + = (1-4) 



=0 for n^d2+q2 
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where the last term vanishes for n ^ q2 +d2 because the bilinear concomitant kernel B2{x; y, s) is a polynomial 
in y of degree d2 — 1 and the linear combination tpn contains the orthogonal function ij^n-qi- 
For the differential equation we have (by definition of the V'n's) 

l^n := V'n + (1-6) 

1 

We then have 

dM' = I JJ e— (5. - ,) g2(x^2^^) eC.-v.fa)+v.(.)^„(g) = 

y - s a; 

^2(y) -^2(s) ^g„_^,_y,fa^+y,(,-) V^n(0 , 

y-s a; 

_ 5)e«''-^«-^=(^)+^=^(«)V^„(0 = 

X — 4 



■2^,(n+j>i7^.+ 



-1 

^ C CC g2(y) -g2(g) ^g^_^s-y.faHV2(s) ^»(0 I 

Jr JJ>t ^ ~ ^ 

di 



-1 



+ Jl (a;(B2(2/) -S2(s)) -^2(2;) +5^(5) +A2(s))e««-^^-^=^('')+^=(^)-'^"^^^ 



(the s-part is a total derivative) 



■2^,(n+j)^S^.+ Ml (xi32(y)-A2(2;))e« 



-1 

di 



r JJ>t 



L 11 S2(y)e«^-^^-^^(^'+^^(^)^n(0 + 



total derivative in y 



{^B2{y) - A2{y)y 



t,P^{n+3)i^'ll, + f ff S2(y)e«^-^^-^^(^)+^^(^)^„(0 (1-7) 



=0 for n^92 + l 

In the step marked with * we have performed an integration by parts: in this integration we do not get any 
boundary contributions because the quasipolynomials ipn by definition are divisible by Bi (which vanishes at 
all endpoints and in particular at the hard-edge ones). This concludes the proof. Q.E.D. 
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A. 2 Proof of Theorem 4.2 

During this and following proofs we use the notation 

^n(y) := , (1-8) 

for the row-vector of quasipolynomials in y. Moreover, at the risk of marginal confusion, we omit all differentials 
of the integration variables since which variables are integrated and on which contour should be always uniquely 
determined by the context (the formulas become significantly longer otherwise). For eq. (4.1) we have (recall 
that A(^) is linear in ^) 

n-l 



(LHS of 4-28) = e""''"^ \\ 0"^^ ^^^^ (x) (z - x) = 



=6 



JJ e-^i«)+«^$„(y)A(e)$„(^) 



1 



AeVi(.) rr e-^UO+«?.£!iMA(;^)$„(a:)-e^^(^) ff e-^^(«)+«^$„(2/)A(x)H^ = 
=$i°)(z)A(z)§„(a;)+e^^W 2 e-^^(«)+«'y0,(2;)^,(x) = 

=ti°'(z)A(z)*„(a;) -Fe^i("'-^i(^' , (1-9) 



where in the identity marked * we have used the linearity of A which implies the following identity 

A(0 A(0 _ Hz) A{x) 
z-S. x-S, z-^ x-£, ' 



(1-10) 



The second form of (a) is proved along the same lines using the principal CDI for the kernel Kn (in Thm. 
4.1). For the remaining GDI's we have 



(LHS of 4-29)=^ "2 e^V.(y)J^ JJ B2{x;^,s} 



Qrip-xs + V2is)-V2M 



V'i(p) 



X — p 

~f e^'^n{y){ ff g.(x;^,.)e''^— .^^"^("~"\ $„(p) = 
2*71- Jr^. JfJJ^ {z-pKx-p) 

— f e--^$„(2/) r rr 62(x;77,s)e'"'-^^+^-(^)-^^('')A(z) ^) $„(p) = 

=$W(z)A(z)$«(x)-^ 2 e^>,(y)J^ JJ S2(:r; 7y, 5)0""— +^^(^)-^^ ^^.(p) = 

=t</)(z)A(z)$„(x)--^ r f S2(x;j/,s)e^^--^+^^(^)-^=(^) , (1-11) 
2«7r Jr, Jf, 
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where the identity marked * is valid for n ^ d2 (so that the kernel reproduces the polynomial 62(2;; ?7,s) of 
degree d2 — 1). 

The proof of the second form of (b) is only marginally different in that we have to use the second form of 
the principal CDI for the kernel Kn (in Thm. 4.1). Q.E.D. 

A. 3 Proof of Theorem 4.3 

Let n — 1 ^ m ^ n + d2 — I: in the following chain of equalities all the steps are "elementary" and hence the 
computation is straightforward. For reader's convenience we have tried to make annotations on the formula 
in order to highlight less obvious steps. 



^^y,(.) rr (l-/(x)-y/(0)e-^^(^)+^V„(;;) ^ ^y^^,, ff y e-^^ffl+^^^</>„ 



(y) 



:(C) 



{-B + C)+ J] 0f (x) ^,(p)7?0„(,7)e 



lb — _L n-\-d2 pf 

j=0 j=n 
n+d2 n — 1 

-i-B + C) + 2 ^^i^)P,m + 2 ifHx) [P,(p)^-\p) 

n — 1 



n+d2 



i-B + C)+ 2 (^)i^.rn + 

j = n 

n — 1 pp 

+ J]^fHx) 0,„(r/)c'"'-^^('')y/(p)^,(p) 

n-\-d2 

(C)+ J] 0\x)P,^ + 

JJ d^.nivV'^Vlip) 



$i")(a;)A(a;)$„(p) +e^i(^)-^iW 



X — p 



+ 



+ 



$1°)(x)A(x)§4p) 

X - p -™ 

V „^^i"' (x) A(rg)§„(p) +c^i( ")-^i (") 
X — p 
^^\x)A{x)^,,{p) 

X - p -™ 

,$W(a:)A(a:)$„(p) 



X — p 
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n+d2 



n-\-d2 



X — p 

$i°'(a;)A(a;)$„(p) + c^i(^)-^i('') 
X — p 



fin 



n — 1 pp 



vlip) - vijx) 

p — X 



n+d2 



2 ^f\x)Pj^ + 



X — p 



ip) 



+ V({x)^^l^{x)S„^,n-l 



$i0Hx)A(x)JJ $„(p)0™(77)c 



p — X 



(1-12) 



We note that in this last expression we have dx(f>^^''{x) expressed purely in terms of (f>^^\x) for £ 



1, . . . n + (^2, the value £ = n + d2 entering only in the first expression. Given that (j)^-^^ (x) satisfies the same 



multiplicative recurrence relations as the Fourier-Laplace transforms for n ^ 1, we can re-express 0^"^^ m 
terms of the elements of the window ^^['^ (a:), obtaining the result. 

The computation for the Fourier-Laplace transforms gives also the same differential equation, indeed 



p n + d2 r^r* 

'n,-\-d-2 n — 1 

2 l^;\x)p,^ + J]l^;Hx) e'"'-^^W0™(,7)(-5p + n'(p)V.(p)+ 

n-1 



j=n 
n-1 



+ 

3=0 
n+d2 



2 </.f)(x)P„„+$W(x)A(x) c"''^™(r;)^vl/„(p) + 



+ ^l:\x)Aix 

n+d2 



^n{p)ct>^:\p) 



X — p 



2 ^f)(2^)^.™+ti'''(^)A(x)JJ $„(p)c'"'0™(r7)^ 



n+d2 



Vi{x)^:\x)K{x) W l!il^e'"'0„(r,)+$(:){x)A(x) 

■'-rd2 pr. 

^f^(x)F,™+$M(x)A(x)JJ $„(p)e''''</'™(77) 



X — p 

^n{p)^::\p) 

X — p 



PS< j;>f 



■IL — 1 pp 

+ n'(^) 2 '(^) ^■,(p)e'"'0™(77) +$('■) (x)A(x) 



X — p 

' ^n{p)^:^Hp) 

X — p 
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Vl(p) - V{ix) 



X — p 



+ Vlix)S^,„-il^\x) +^^\x)Aix) 

The coefficients of these expressions in terms of 
computation, hence completing the proof. Q.E.D. 

A. 4 Proof of Theorem 4.4 

Let n — d2 ^ m ^ n and let us compute 



X — p 



n — 1 —n-\-d2 — l 



(1-13) 

are precisely the same as for the previous 



n—r\ JJ>c 
n — 1 



n-1 

V{{x)-4)jn{x) + 2 -Ipjix 



Vi{x)4!m{x) + 
Vi{x)4>,n{x) + 



i-x 



n-1 



n-1 



i = m-l -^-^^ 

VY(2^)^mnV'n(2^) + 



IL — J. 

A(x)§,(x) - 2 ^.(2^) - Vlip))$,iy) 

j=0 
n — 1 

A(x)^,{x) + 2 ^j(^) i^miO<^^"Vlip)$,(y) + 



Hx)^,{x) + 



n — L pp n — i 

j=0 ^-^^ j=m-l 



Vl{x)6mni'n{x) + 



^-X 



A{x)^,{x) + 



.iu n{p)-v({x) 



p — X 



$n(2/)A(x)§j(x)- 2 ^j(^)^" 

j = m-l 



(1-14) 



The last term contains ^n-d2-i (for m = n — ^2) which is "outside" of the window of the quasipolynomials. 
Using the recurrence relations and re-expressing it in terms of elements in the window (using the ladder 
matrices) we obtain the formula. 

For completeness one should also consider the other columns of the fundamental system '^n and show 
that they satisfy the same differential relation as the quasipolynomials. Let n — d2 ^ m ^ n, then 



5. 



V'm(e) 

x-C 
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x-i 

B2iy)-B2is) B2{x;y,s) 

{y-s){x-0 ' {x-0 



B2ix;y,s)8^ 



x-^_ 



= :(S) 



,iy-xs + V2is)-V2{y) V'm{0 
iy-xs + V2{s)-V2{y) 



x-i 

B2{y)-B2{s) 



2/ - s 

,Cy-2;s + V2(s)-y2(y) (55-v\'(0)p™(0 
x-S, 

hoSmO because V2— V2— lnS2 



+ (2/ - s)B2{x;y,s) 



2«7r Jf „ JJ^ a; - ^ 



+ B2{x;y,s)i 



,Cy-2;s+V2(s)-y2(y 



:(Z5) 



(B)-(C)+VM: 



-Vi(i)+iy 



{B)-{c) + {D) + Y,rp{ 



j=0 



j=0 



n-1 



Ux){l-5ra.,nW^Kx) - i^f\x) i^^OyUv)^ 

-1— n JJ>c 

A(x)§„(x) + 



+ 1^1 



j=0 



i-x 



^-X 



B2{x;y,s)e 



yi-xs + V2{s)-V2{y) 



[gives (C)] 
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IL — J. 

j = m-l 



j = m-l 



n-1 



<5„„F/(x)^M(a;) - ^ \\ i^miOvMy)^^' (1-15) 



3=0 



This is the same expression as for the quasipolynomials: since the auxiliary wave functions 'ipj^\x) satisfy the 
same multiplicative recurrence relation (for n large enough) as the quasipolynomials, re-expressing ipnliix) 
in terms of the elements of the window yields the same differential equation. Q.E.D. 

A. 5 Proof of Theorem 5.1 

For brevity we denote A„(x) simply by A{x) during this proof. Since the rows (columns) of 4>„ (^n) are of 
two types, we need to carry out four types of computations 

(a) = $W(x)A(x)*l°'(x) , (b) = $W(a:)A(a:)$«(x) , j = l...d, 

(c) = {x)A{x)¥°^ (x) , J = l...d2 (d) = $ W (x) A(x)$(r^ {x) , £,m = l...d2 (1-16) 

It follows trivially from (5-3) that (c) = (set x = x' in the LHS). For (a) we have 

(a) = e^^(-) \[ i!il§e-^^(«+«A(a:)$„(a;) = (1-17) 

n — 1 



3=0 



where we have used that 4>j{Q, i ^ 1 are orthogonal to p{Q = 1. Note also that we had to use the CDI in 
the form (5-2). Then we have to compute for 1 ^ £,77i ^ d2 (we suppress explicit reference to the variables 
of integration because there is no possibility of ambiguity) 

(d) = rr r S2(a;;77,s)^^^%^e«''--^-^^W+^^(^)= (1-18) 

= ^g^fn^)JJ S2(x;ry,s)V',(0c«''-^^-^=('')+^^(^)= (1-19) 

n — 1 

= 2^2 Jpd.sJ^dy0,(2/)c^^^JJ B2{x;rj,s)i;,{0^i^-^'-''^(^^+^^('^^ ^ (1-20) 



3=0 -^r™ Jri 

— r ds rdy62(x;y,s)e^(''-^)-^^(«)+^^('') = (1-21) 

I32{^y,l,^y,m) = hm , (1-22) 
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where in the step marked with a star we have used that for the polynomial of 77 P{ri) := 62(2;; 77,5) is 
reproduced by the kernel 

n— 1 ™« 

P{y) = 2 ^Ay) dr7dC7^,(Oe-^=W+«''P(r;) (1-23) 

j=Q JJ'^ 

provided that n — \ ^ dcgP = d2 — 1. Note also that in this latter computation we are forced to use the 
other form of the CDI (5-3). Finally we need to compute (b), which involves quintuple integrals 



(b) ^ ^ rr M4e-^(c)+pc rr r ^,(.;,,.)^M^e«----^(^)--^(^) ^ 

- e-^'^^. rr rr r ^,(,;,,.)^(oM)e«.-^^-v.(.).v.(.) + (i_24) 

JJ>. a:^ - C JJ>. Jsef ,„ a:^ - ? 



,{i(-a;s-y2(i;)-l-V2(s) 



=0 if n^qi 



2» |JJlL^"""""'"'""""''''*"'''*'<^*<«''"<-'-'v 



-e' 



rrrds B2ix; r?, 5)^e«"--^-^^ 
2«7r JJJf x-^ 

^ Idy rrdse-^^(«)-^=^('')+';''— +^=^(^'S2(x;p,s)^ + 
2«7r J JJf a:^ - C 



JJJdse-^i(«l-^^('')+«''-"^+^^(^'fi2(a;;?7,s)^^ = (1-25) 



2iTT 

Once more, we are forced to use the CDI in the form (5-3). This concludes the proof. Q.E.D. 
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